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Abstract 

We investigate the scaling properties in inclusive hadron production and the associated single 
transverse spin asymmetry in the forward rapidity region at RHIC. We find that the spin-averaged 
I experimental data in both pp and dAu collisions demonstrates a transverse-momentum-dependent 

■ geometric scaling. We introduce the transverse momentum dependent Collins fragmentation func- 

^ ■ tion to study the scaling of the single transverse spin asymmetries. The general feature of the 

scaling analysis is consistent with the experimental observations, in particular, for the transverse 
I momentum dependence of the spin asymmetries at RHIC energy. We further propose to probe the 

saturation scale of nucleus by measuring the spin asymmetry normalized by that in pp scattering 
'■q^' at low transverse momentum. 
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I. INTRODUCTION 



Since its operation, the forward rapidity region in pp and dAu collisions at the Relativistic 
Heavy Ion Collider (RHIC) has provided important opportunities to study novel hadronic 
physics phenomena, from the single transverse spin asymmetry (SSA) in pp collisions 
to the small- a; gluon saturation in dAu collisions |5|-l8|. Both physics has attracted great 
attentions in the last few years. 

The nuclear suppression observed in hadron production in the forward rapidity region in 
dAu collisions at RHIC has been interpreted as the consequence of the saturation physics 
at small-x [1-14|. An effective theory of the color-glass-condensate (CGC) has been applied 



to describe these phenomena [15[. Alternative approaches have also been proposed |16Hl9 



Meanwhile the large single transverse spin asymmetries in pp collisions found at RHIC have 



been studied from various approaches ^20|, |21j, while the underlying mechanism remains 



unclear [2J]. In this paper, we will extend the CGC formalism to include the spin effects, 
in particular, in the fragmentation process. Early attempts have been reported in Ref . 23| , 
where they focused on the spin dependent quark distribution contribution. 

One important aspect in the context of the color-glass-condensate is the so-called geo- 



metric scaling [2J, |25| . This is better understood as that the unintegrated gluon distribution 
NF{x,q±) can be written as a single function of q\/Ql{x) with Qs{x) the saturation scale. 
As a result, the structure function of deep inelastic scattering can be expressed as a single 
function of r = Q^/Q^ with Q the virtuality of the exchanged photon |25|. It is remarkable 
that the experimental data from HERA support this geometric scaling. Since then, this 
idea has been applied to many hadronic processes involving small- a; parton distributions, 
including some experimental data relevant to the forward region at RHIC. 

In the forward rapidity region of pp (or pA) collisions at high energy, the incoming parton 
from the projectile scatters on the target. Because of high density of the partons from the 
target at small x (or in large nucleus), this parton goes through multiple scattering before it 
fragments into the final state hadron. The CGC formalism takes into account these multiple 
interaction effects in a systematic manner. In particular, the quantum evolution at small- a; 
in this formalism has been well established and plays a very important role to understand 
the suppression phenomena observed at RHIC |26| . 

We further introduce the transverse momentum dependence in the fragmentation func- 
tions to understand the single transverse spin asymmetries by introducing the spin-dependent 



fragmentation function - the Collins function j27|. Previous model studies have claimed that 



this effect is too small to make any significant contribution to the SSA in forward hadron 



production in pp collisions [28[. However, we will find that the CGC formalism provides 
a simple factorized form for the spin- dependent cross section, from which we find that the 
Collins mechanism indeed plays an important role. Of course, the limitation of the current 
knowledge of the quark transversity distribution and the Collins fragmentation function 
do not allow us to make reliable predictions. However, from the positivity bounds of both 
functions, we find that it is possible to explain the large single spin asymmetries in hadron 
production observed by the experiments at RHIC. 

The rest of this paper is organized as follows. In Sec. II, we analyze the inclusive hadron 
production in the forward rapidity region in pp and pA collisions following the CGC ap- 
proach. We will show that the experimental data from RHIC demonstrates a transverse 
momentum dependent geometric scaling for both pp and dAu collisions. In Sec. Ill, we in- 
troduce the transverse momentum dependence in the fragmentation process to understand 
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the single spin asymmetries in pp collisions. In particular, the transverse momentum de- 
pendence of A]sf indicates the importance of details in the hadronization. We conclude our 
paper in Sec. IV. 



II. TRANSVERSE MOMENTUM DEPENDENT GEOMETRIC SCALING 

In the forward rapidity region, the hadron production in pp and pA collisions at high 
energy is dominated by the valence quark and gluon from the projectile. The energetic parton 
penetrates in the nucleon (nucleus) target with multiple scattering and then fragments into 
the final state hadron. We follow the CGC approach to study the differential cross section 
of hadron production p + p(or A) ^ h + X in the forward rapidity region 

da K dz 



dynd^Phi. (27r)2 



f dz 

/ —Xiq{xi)NF{x2,k^ = Ph^lz)Dh/q{z), (1) 

Jxp ^ 

where K is a, possible /^-factor to absorb the higher order corrections, xp = Ph±/ ^/s^^viyh), 
and Xi = xp/z is the momentum fraction of the projectile carried by the incoming quark, 
z the momentum fraction of the quark carried by the final state hadron, X2 = Xi exp(— 2j//j) 
is the momentum fraction of the target participating in the hard scattering. We have also 
suppressed the scale dependence in various factors in the above formula, for which the precise 
forms will depend on the next-to-leading order perturbative corrections. 

Nf{x2, k±) is the unintegrated gluon distribution, and is given by the two-dimensional 
Fourier transform of the imaginary part of the forward dipole-target scattering amplitude 
in the fundamental representation: 

Nf{x2, k^) = j dhe-'"^-' [1 - A/>(r, Y = \n{xo/x))] , (2) 

with r the dipole size. The characteristic feature of A/'F(r, Y) is its "geometric scaling", i.e.. 



J\fF{r,Y) only depends on the dimensionless quantity r'^Q^{x) |2J, |25|. Thus its Fourier 
transform has "geometric scaling": k'j_NF{x2, k±) depends only on the dimensionless ra- 
tio kj_/Ql{x). This property will manifest itself in the differential cross section of hadron 
production as we will show below. 

For the large- a; quark and gluon distribution, an important power counting can be app lied 



to study the general behavior for the differential cross section of the above process [29|, |30 |. 
In particular, for the region where the quark distribution dominates, we have q{x) ~ {1 — x)^. 
Furthermore, fragmentation function also has power counting rule, Dh/q{z) ~ {1 — z), where 
the power behavior comes from the current parameterization of the fragmentation function 
for the charged and neutral pions Combining the above power counting analysis, we 

obtain the following power behavior for the differential cross section, 

^-dk - A)) ■ 

where the saturation scale Qs depends on the target and the momentum fraction X2 carried 
by the gluon. The geometric scaling (as a function of Ph±/ only) comes from the geometric 
scaling of the unintegrated gluon distribution for k\NF{ mentioned above. 

^ The original power counting would predict (1 — z)^ behavior for the spin-averag fragmentation func- 
tions [30|. This would lead to a power behavior of (1 — xp)^ instead of (1 — xf)^ in Eq. ([3]). The current 
data in Fig. 1 can be described by both parameterizations. Future experimental data at very forward 
region shall be able to distinguish them. 
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Similar geometric scaling for inclusive hadron production at central rapidity has been 



observed in [32|, |33|], where it is the differential cross section alone exhibits the geometric 



scaling. On the contrary, we have to multiply a factor of , to the differential cross section 



to show the geometric scaling. The geometric scaling derived in [32|, |33| relies on the fcy- 
factorization (implicitly), and both incoming hadrons contribute to the saturation effects in 
the particle production at mid-rapidity, where the saturation scales from the two incoming 
hadrons are in the same order at mid-rapidity. However, in our analysis, we have a dilute 
projectile scattering on a dense target, where the saturation only comes from the dense 
target and the geometric scaling is manifest in the first place. Of course, as a result, the 
differential cross section has to be multiplied by a factor of P^^ to demonstrate the geometric 
scaling. 
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FIG. 1. Geometric scaling of the differential cross sections (left) and the nuclear modification 
factor i?dAu (right) plotted as a function of P^^/Q^. Data are from BRAHMS and STAR " 
Collaborations at RHIC. 



In Fig. [H we plot the RHIC data in pp collisions the differential cross section, 

2 da{pp h + X) 



[i-XFr'p, 



dyh(PP, 



hi. 



(4) 



as a function of Pl^/Ql{x2 



the saturation scale Ql{x2) 
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where X2 is calculated as X2 = P^^^/^Xis) ~ P^j_/ (xpS), and 
- Ql{xo/x2)^ with Qo = 1 GeV, Xq = 3 ■ 10-^ and A = 0.28 
In this plot, we have included the data for vr^ from BRAHMS collaboration jsf and 
from the STAR collaboration ^ at y/s = 200 GeV. tt^ data are from two rapidity bins 
Uh = 2.95 and i/h = 3.3, whereas those for 7r° cover i/h = 3.3, 3.8,4.0. We have also included 
data for 7i~ production at lower energy ^/s = 62.4 GeV from BRAHMS at i/h = 2.7 and 
yh = 3.3 jsl. For 7r° production, we have multiplied a factor 2.5, which is consistent with 
the overall i^-factor K ~ 0.4 found in 13||. From this plot, we see clearly that the hadron 
production in the forward rapidity region demonstrates the geometric scaling. 

It is also interesting to check the geometric scaling when we compare the pA and pp 
collisions in the forward rapidity region 



34i |. Here, the nuclear modification factor RpA is 



measured, 



N,oiidNPP^f^+^/dyhd^Phi. 



(5) 
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where Ncoii represents the number of binary colhsions in the selected nucleon-nucleus scat- 
tering. Many factors cancel out in the above ratio. Therefore, according to Eq. ([3]) we will 
have the following geometric scaling behavior, 

'''' 

where again the saturation scale Qs depends on the momentum fraction of the nucleus 
target involved in the scattering process and is determined by the kinematics of Ph± and 
rapidity i/h- In the above equation the subscript b indicates that the ratio also depends on 
the centrality of the collision, which determines the change of the saturation scale from pp 
collision to pA collision at a given impact parameter. Of course, the magnitude TZ depends 
on the overall normalization as well, including the parameter Ncou- 

As an illustration, we plot the -RdAu for hadron production at ^/s = 200 GeV in the 
forward rapidity region, where we include the BRAHMS data for the negative charged 
hadron h~ at two rapidities: yh = 2.2 and i/h = 3.3, and the STAR data for neutral pion at 
Hh = 4.0. Clearly, the two rapidity sets from BRAHMS demonstrate the geometric scaling 
behavior, whereas the STAR data are a little off the scaling curve of the negative charged 
hadrons. This indicates additional mechanism for very forward hadron suppression in dA 



collisions [19|. It will be very interesting to further test this geometric scaling in the future 
experiments at RHIC and the LHC 

The power behavior for the differential cross section has been observed in various scatter- 
ing energies, and has been used as a practice to understand the forward hadron production. 
However, the geometric scaling only appears when we hit the saturation region at small-x 
of the target. For example, we notice that the forward hadron production at low energy 
fixed-target experiments does not show the same scaling behavior jsHj. This tells us that 
geometric scaling is not yet reached in these energies, and the production mechanism differs. 
Similar conclusions have been drawn from the comparison between the RHIC data and those 
from the previous fixed-target experiments in the coUinear factorization approach j36|, 37 



III. SINGLE SPIN ASYMMETRY IN THE FORWARD REGION 

In this section, we extend the above formalism to calculate the single spin asymmetry 
in the forward region. For this purpose, we need to introduce the transverse momentum 
dependence either in the incoming parton distribution from the projectile (polarized pro- 
ton) j23j and/or the fragmentation function. Here we concentrate on the contribution from 
the fragmentation function. 

In the forward rapidity, introducing the transverse momentum dependence in the frag- 
mentation function in Eq. ([T]), we have the spin-averaged differential cross section 

^ [ -J [ d'^PhTXiq{xi)NFix2,kj_)Dh/qiz,PhT), (7) 



dyhd^Ph± (27r)2 



where PhT is the transverse momentum of the final state hadron relative to the fragmenting 
quark jet, thus the momentum of the final state hadron in Lab frame can be written as: 
Ph± ^ zk± + PhT- On the other hand, the spin- dependent differential cross section for 



p^ + p(or A) ^ h + X can be written as [38 
dAa K dz 



dyhd^Pk± (27r)2 / z^ 



d PhtHSj., PhT)xih{xi)NF{x2, k±)6q{z, Put) 
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where h{xi) is the quark transversity, Sq{z, Phr) is the Colhns function that is related to the 
Trento convention j39| as: 5q{z,PhT) = —H^/zMh with Mh the final-state hadron mass. 
The azimuthal angle dependence is encoded in /(S"^, Pht) and is given by 



I{S^,PHT) = e^^Sl 



p/3 _ ^-PhT p/3 

n-Pj ' 



(9) 



where e"^ is the 2-dimensional Levi-Civita tensor with e^^ = 1, n is a unit vector along the 
momentum of the unpolarized hadron, and Pj is the momentum of the fragmenting quark 
jet. In the forward rapidity region, I{S±, Phr) can be reduced to the following form, 

I{S^,PhT) = |^±||P/.T|sin(0;,-0,), (10) 

where 0s is the azimuthal angle for the spin vector S± in the Lab frame, (ph is the azimuthal 
angle of the produced hadron momentum Ph± relative to the fragmenting quark jet. 
The single transverse spin asymmetry is defined as 

_ dAa / da , , 

dynd^Ph^/ dynd^Ph^' ^ ' 

To be consistent with the experimental definition we choose the coordinate system where 
hadron transverse momentum Ph± along x-direction, the spin vector S±_ along y-direction, 
and the polarized beam along z-direction. 

Before we present the numeric estimates, we would like to study the generic feature 
for the single spin asymmetry in the above approach. We will take a simple model for 
the unintegrated gluon distribution and the Collins fragmentation function to study the 
transverse momentum behavior of the single spin asymmetry. First, we will examine the 
single spin asymmetry at relative low transverse momentum at order of the saturation scale: 
Phi. ~ Q% We know that the unintegrated gluon distribution has saturation behavior at 
low transverse momentum. As a convenient model, we follow the GBW parameterization 
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iVF(x,g^)~-^e-''i/«^ (12) 

which has geometric scaling. For the fragmentation function we choose a simple Gaussian 
distribution, 

Dn,,{z,p^) ^ ^^e-^y-^ (13) 

where represents the width of the transverse momentum dependence in the fragmentation 
function. By performing the transverse momentum integration in Eq. ([7j), we will find that 
the spin-averaged cross section has the following transverse momentum behavior, 

1 -i^?^/(Q2+A2)^ ^^^^ 



dyhd^Ph± Ql + A 

where we limit the above results at small transverse momentum region of order of Qs- If the 
saturation scale Qs is much larger than the transverse momentum width A in the fragmenta- 
tion function, we will find that the transverse momentum dependence in the fragmentation 
function does not change the geometric scaling discussed in the last section. 
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Now we turn to the spin-dependent differential cross section. For the Colhns function, 
we assume it also has a Gaussian form (38| . 

with a slight difference in the Gaussian width — 5"^ to satisfy the positive bound. Substi- 
tuting the above equation to the differential cross section formula in Eq. ([8]), one obtains, 

From the above results, we find that the single spin asymmetry behaves as, 

I A 



,2 -e"^ , (17) 



p. , A ^''■Pfex 



where we have made reasonable assumptions: Q^g ^ A^ 5^. The above result indicates 
that the asymmetry vanishes when P^x — 0, and it also depends on the transverse momen- 
tum width in the fragmentation function. Certainly, if there is no transverse momentum 
dependence, the whole effects will vanish. Furthermore, the spin asymmetry also decreases 
with the saturation scale. This is because the fragmentation effect is suppressed if we in- 
crease the saturation scale (see also Eq. (fT^ ). From the above simple analysis, we find that 
the spin asymmetry in general will have broader distribution as a function of Ph±- 

Moreover, it is interesting to note that the double ratio of the spin asymmetries between 
A and p^p scatterings scales as 



AP 



APP- 



2 d2 x2 
2 



(18) 



at small transverse momentum, where we have assumed that the saturation scale for nucleus 
is much larger than that for the nucleon at the same kinematics. This is the most interesting 
result from the scaling analysis. The ratio of the spin asymmetry is inversely proportional 
to the saturation scale in the limit of Ph± — )■ 0. This can be used as an important probe for 
the saturation scale of the gluon distribution in the nuclei target. 

Similarly, we can estimate the large transverse momentum behavior for the spin asym- 
metries, where the unintegrated gluon distribution behaves as 

iVF(x,gx)~^. (19) 
(11 

If we still assume that the fragmentation function can be parametrized as a Gaussian form, 
we will find out, 



AN{Ph±) OC ^ , (20) 
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where the factor 4 comes from the power of the unintegrated gluon distribution at large 
transverse momentum. The asymmetry decreases as 1/Phi. at large transverse momentum 
as expected. However, the rate of the decreasing is strongly affected by the relative size 
between P^^ and A^. This additional 4A^ modification is slightly different from the usual 
1/Pft_i_ power counting results [20]. It comes from the effects of the transverse momentum 
dependent fragmentation function in the spin-averaged cross section which was neglected 
previously (20| . 

Furthermore, we notice that the saturation scale dependence cancels out between the 
spin-averaged and the spin-dependent cross sections, and the asymmetry does not depend 
on the saturation scale. As a consequence, the double ratio will approach 



A- 



pA—i-h 



APP- 



^ 1 . (21) 



The above results are very interesting observations, and it will be important to test these 
predictions in the future experiments (4l| . 




FIG. 2. Left: Ph i -dependence of the single spin asymmetry. We use Eq. (fT7|) to fit the STAR 
experimental data [^jli^]- Right: We plot Ajye^^ ^h±/Qs Q(^+'^^)y g, function oi xp - the generalized 
j;F-scaling for the single spin asymmetry according to Eq. (j23p . We have chosen J = 0.16 GeV to 
demonstrate the scaling. The data is from STAR ^]. 

Comparing Eqs. ( IT7|) and (!20|) . we find that the maximum of the single spin asymmetry 
will be at P^^ ~ Q^/S"^. Since 5^ is a small number, we would find that the spin asymmetry 
drops at a relatively large transverse momentum, which seems to be consistent with the 
observations at RHIC experiments 0, E^. In Fig. m (left), we show the comparison of 
Eq.( IT711 with the experimental data from the STAR collaboration [2I, |42[, where we have 
chosen ^ 1.0 GeV^, which is roughly the average saturation scale in this experimental 
kinematic region. 6 ~ 0.16 GeV is adjusted to fit the experimental data. Although the 
normalization of the asymmetry depends on and the size of the Collins effect, including 
the quark transversity distribution and the Collins fragmentation function, the transverse 
momentum dependence from our general analysis is consistent with the experimental data. 

Meanwhile, from Eq. f|T7|) and realizing Ph± = XFy/Se"^'^, we find 
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where T [xp) represents any additional xp dependence for the single spin asymmetry. There- 
fore, the following quantity 

A;ve^'^-'-/«^e(^+2^)^'' ~ x^^^'^T (xf) , (23) 



will be a universal function of xp, independent of the rapidities. In Fig. [2] (right), we plot 
STAR 7r° single spin asymmetry data for the quantity defined in Eq. (123!) as a function of 
xp for two rapidities yt = 3-3 and yt = 3.7. We find that for relatively small S ^ 0.2 
GeV (which is consistent with our approximation), the data seem to be consistent with the 
Xi7-scaling. It will be interesting to test this scaling in the future experiments. Note that we 
derive this interesting scaling purely from the fragmentation process, but it could also come 
from the distribution function (Sivers function) in the polarized nucleon However, the 
complicated structure for the Sivers function makes such a study nontrivial (22| . 




FIG. 3. The comparison with the STAR experimental data at ^/s = 200 GeV for rapidity yh = 3.3 
(left) and yh = 3.7 fright). The solid lines are using the transversity from 43|, |4J] and the Collins 
function from 38|, |4J]. The dashed curves are the upper bound as explained in the text. 
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FIG. 4. The comparison with BRAHMS experimental data at ^/s = 200 GeV. The curves are the 
same as explained in Fig. [3l 
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Let us now come back to make some numerical estimates for the single spin asymme- 
try. For this purpose, we will use the unintegrated gluon distribution in Ref. [l3(, which 
comes from solving B al it sk y- Kovchegov equation with running-coupling corrections. We will 



choose the quark transversity in |43l . |44| , and adopt the transverse momentum dependent 



unpolarized fragmentation function and Collins function in jssl . Si 



Dh/q{z,PhT) 



6q{z, Phx) 



n{pl) 
2 



h/q{Z), 



(vr(pi)) 



2 \n3/2 



e-nV/(pi)5g(V2)( 



(24) 
(25) 



where (pj) = 0.2 GeV^, and 5q^^^'^\z) is the so-called half-moment of the Collins function 



given m 



The numerical estimates for the single spin asymmetry of vr^ are shown in 
Fig. |3]as solid lines, and also compared with the experimental data from STAR collaboration 
[2|. Similar comparison for the charged pions are shown in Fig. Hlwith the experimental data 
from BRAHMS collaboration fl. 

From these plots, it seems that the predictions are below the experimental data. However, 
it is important to keep in mind that there are uncertainties in both the quark transversity 
distribution and the Collins fragmentation function. For example, the quark transversity 
is constrained using experimental data from HERMES collaboration, only up to s ^ 0.3 



45[, thus we are lacking of information on the transversity in the interested region at RHIC 
forward direction x > 0.3. The situation for the Collins function is slightly better since 
we have complimentary information from both semi-inclusive deep inelastic scattering at 



HERMES |45| and e'^e~ collisions at BELLE |46|. But one still needs to understand how 
to incorporate the evolution of the Collins function, which will affect the Gaussian approx- 
imation and the Gaussian width used above [47|. Realizing the limitation of the current 
knowledge for both quark transversity and Collins function, we will set up a bound for the 
single spin asymmetry by using the positivity bounds for both functions. 

For transversity, we saturate the Soffer bound |4^], which is derived from positivity j49| . 



Following Ref. j28|], we choose 



1 1 

hu{x) = - [u{x) + Au{x)] , hd{x) = -- [d{x) + Ad{x)] , 



(26) 



where Aq{x) is the helicity distribution function, hq{x) is the quark transversity. For the 
Collins function, we have 



PHTSd^\z,PHT)=Df{z,PhT). 



(27) 
(28) 



where the signs are chosen such that we have the maximum single spin asymmetry for the 
charged pions [28!]. For 7r°, we rely on the isospin symmetry. In Figs. [3] and HI the dashed 
lines are the upper bound for the single spin asymmetries. From these plots, we see the 
Collins mechanism in CGC formalism is large enough to explain the observed asymmetries. 



If the Collins mechanism is the dominant contribution to the single spin asymmetries |22 



one should be able to extract the quark transversity and Collins function by using the RHIC 
data. 
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IV. SUMMARY 



We study the scaling properties in the inclusive single hadron production and the asso- 
ciated single transverse spin asymmetries in the forward rapidity region at RHIC. For the 
spin-averaged differential cross sections in both pp and pA collisions, we find a transverse 
momentum dependent geometric scaling. These scalings are clearly seen in the spin-averaged 
experimental data in both pp and dAu collisions. We further introduce the transverse mo- 
mentum dependence in the fragmentation function to investigate the scaling of the single 
transverse spin asymmetry associated with the Collins mechanism. The general behavior of 
the transverse momentum dependence of the single spin asymmetry is consistent with the 
trend observed in the experiments at RHIC. We also find that the double ratio of the spin 
asymmetries in pA and pp collisions is inversely proportional to the saturation scale in the 
limit of Ph± — 0. This can be used to probe the saturation scale in the nucleus by measuring 
the spin asymmetry normalized by that in pp scattering at low transverse momentum. 

With the limited knowledge for the relevant functions - quark transversity and Collins 
function, we set up an upper bound for the single spin asymmetry. We find that the bound 
is certainly large enough to explain the experimental data at RHIC. 
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